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duced in order to describe the behavior of neutrinos in matter. The former 

relates weak eigenstates to mass eigenstates, whereas the latter relates weak 

eigenstates to energy eigenstates in matter. It is shown that the dynami- 
cs ■ 

S-h ' cal mixing angles enable us to express the neutrino survival probability in 

the Sun without any resort to the Landau-Zener transition probability for 
the non-adiabatic process. Also discussed are effective CP violating phases 
that appear in the effective and dynamical mixing matrices in matter. Both 
two and three generation cases are discussed using the solar neutrinos as an 
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I. INTRODUCTION 



The Mikheyev-Smirnov-Wolfenstein (MSW) effect has been considered for some 
time as the most plausible mechanism to explain the solar neutrino puzzles |§-|6|]. This 
effect dramatically enhances the conversion of a flavor neutrino [y e in the case of solar 
neutrinos) into another when certain conditions among the intrinsic neutrino properties and 
the nature of matter are met in a medium. When the matter density in the Sun changes 
slowly enough (to be determined by the intrinsic properties such as masses and mixing angles 
and by the neutrino energy) so that the neutrino conversion from one flavor to another takes 
place adiabatically, the MSW effects are known to be well described by the effective masses 
and mixing angles in matter. These effective values at any given time or distance from 
the production in matter are commonly derived from the equation of motion for the weak 
eigenstates in matter by a simple matrix diagonalization under the assumption that the 
matter density is uniform. 

However, when the time (or distance) variation of the density in matter is sufficiently 
rapid, transitions among the effective mass eigenstates take place and the process becomes 
non-adiabatic. In this case, the role of the effective values becomes vague. In order to 
overcome this, the neutrino survival probability is approximately expressed in terms of the 
Landau-Zener transition probability |7Hl4|. The other alternative is, of course, to solve nu- 
merically the equation of motion for the weak eigenstates in matter. It has been recognized, 
however, that to do so is prohibitively difficult when all three generations of the neutrinos 
are involved. 

In this paper,we first introduce two concepts of the mixing angles in matter, effective and 
dynamical mixing angles and the corresponding mixing matrices. It will be shown that the 
dynamical mixing angles enable us to calculate the survival probability without resort to the 
Landau-Zener transition probability. The effective mixing angles defined here are equivalent 
to those that have been commonly used in the past, i.e., they appear in the effective mixing 
matrix which relates the weak eigenstates and the mass eigenstates in matter. As mentioned 
already, these angles are sufficient to describe the adiabatic process and have also been used 
for the non-adiabatic case together with the Landau-Zener transition probability. 

The dynamical mixing angles are defined to appear in the mixing matrix which relates 
the weak eigenstates and the energy eigenstates in matter. (Note that even in the case 
when there appear transitions among the mass eigenstates, the energy eigenstates are well 
defined.) It will be shown that the neutrino survival probability can be expressed in terms 
of the dynamical mixing angle alone even in the case of non-adiabatic processes. Since in 



this approach no approximations are involved in calculating the survival probability, the 
results are exact whereas the usual treatment with the use of the Landau-Zener transition 
probability relies on some approximations. 

The second subject to be discussed in this paper is the medium effect on the CP violat- 
ing phases already present in the neutrino mixing matrix. In the case of general mixings, 
it has been known that the CP phases would not significantly modify the neutrino transi- 
tion probabilities, in particular when the oscillating terms are averaged out. Therefore, an 
observation of the CP violating effects is rather difficult, if not impossible. This general 
conclusion, however, needs not to be correct in the case of maximal mixings. The possi- 
ble solution of the atmospheric neutrinos and the so-called large angle solution of the solar 
neutrino problem require rather large mixing angles, suggesting the possibility of maximal 
mixings of three neutrinos • In this case, it is necessary to have a CP violating phase. 

That is, maximal mixings are impossible without CP violating phases jHJ . 



In this paper, we investigate how the medium effects would modify the CP violating 
phases in the mixing matrix in the vacuum. It will be shown that in the case of three 
generations of neutrinos, in addition to numerical changes in the original vacuum CP phase, 
the medium effects induce additional effective CP phases which, of course, would vanish 
when the matter effects are turned off. Although it is necessary to have additional CP 
phases in matter, they are not independent of the original vacuum CP phases. In the two 
generation case, no additional CP phase appears in matter. 

This paper has been organized as follows: In Section II we introduce two mixing matrices, 
the effective mixing matrix and the dynamical mixing matrix. Section III is devoted to a 
detail discussion of the effective mixing matrices in two and three generations of neutrinos. In 
particular, we demonstrate how additional induced phases appear in a medium in the case of 
three generations. In Section IV, we discuss the dynamical mixing matrices. After presenting 
a simple example of the two generation case, we discuss the case of three generations with 
non-zero CP. We summarize our results and conclude in Section V. 



II. DEFINITIONS OF EFFECTIVE AND DYNAMICAL MIXING MATRICES 

Let us begin with the wave equation for neutrinos propagating through a medium with 
a time-dependent potential, V(t), being felt by the neutrinos a la MSW effects 

^||<M^) >= ^M^ w (t) >, (1) 



where \^ w (t) > is the neutrino wave function in the weak (flavor) basis and M w is the 
time evolution matrix in the weak basis which dictates the motion of the weak eigenstate 
neutrinos. First, we define the effective mixing matrix U e ff which satisfies the following 
equation 

Ut ff M w U eff = M"\ (2) 

where M m is the diagonal mass matrix in the mass basis. Solving the neutrino wave equation 
Eq. ([I]) is equivalent to finding the effective mixing matrix U e ff in Eq. (0) in the case where 
neutrinos propagate in a medium with a uniform density (V(t) = constant). The result is 
sufficient to describe the adiabatic MSW process in a medium. 

Next, we introduce a unitary matrix U which diagonalizes the time-evolution matrix for a 
time-dependent V(t). By transforming the wave function in the weak basis with the unitary 
matrix U, we obtain the wave function in the energy basis such that its time evolution is of 
the form 

|* e (t) >= e-^/o M W|# e (0) > (3) 

where M e (t) is diagonal and the subscript e in the wave function \l/ e denotes the energy 
eigenstate. In other words, the above-introduced unitary matrix U relates \^ w (t) > and 
|\l/ e (0) > for a time-dependent V{f) as 

\* w (t)>=U(t)\V e (t)>. (4) 

By inserting Eq. (|j) into Eq. (|l]), we can obtain a differential equation for the matrix U . 
From now on , we will call this U the dynamical mixing matrix and use the notation Ua yn . 
The Udyn satisfies the differential equation 

U\ yn M w U dyn = M e + 2EiUlJj dyn , (5) 

where M e is a diagonal matrix whose elements are composed of time varying energies of the 
Hamiltonian M W /2E. 

The form of the dynamical mixing matrix Ud yn does not depend on whether the neutrinos 
are Dirac or Majorana type. For the two generation case, the dynamical mixing matrix is 
a 2 x 2 matrix with one dynamical mixing angle and one dynamical phase if the effective 
potential V(t) is not a constant in time. For the three generation case, the dynamical mixing 
matrix is a 3 x 3 matrix with three dynamical mixing angles and three dynamical phases. 

Summarizing the above, the effective mixing matrix relates the weak eigenstates and the 
mass eigenstates 



(6) 



whereas the dynamical mixing matrix relates the weak eigenstates and the energy eigenstates 

| weak eigenstates >= Ud yn \ energy eigenstates > . (7) 

Usually, it is easy to calculate the effective mixing matrix by using simple algebra , i.e., 
it is an eigenvalue and eigenvector problem. The eigenvalues correspond to the effective 
masses and an appropriate assemble of eigenvectors is nothing but the effective mixing 
matrix defined by Eq. As for the dynamical mixing matrix, however, it is non-trivial 
to obtain the dynamical mixing angles, phases and the time-dependent energy eigenvalues 
because it is equivalent to solving the Schrodinger equation of the neutrino wave function 
given in Eq. ([!]). Therefore, one often resorts to numerical methods to solve Eq.(l). 



III. EFFECTIVE NEUTRINO MIXING MATRICES 



A. Two Generation Case 



As a simple example, let us consider the two generation mixing of Majorana neutrinos to 
investigate how the CP phase in vacuum would be modified in matter. (In the case of Dirac 
neutrinos, the phase is trivial, that is, zero in matter as well as in vacuum.) In general, the 
mixing matrix in vacuum can be written as 



o 



cos 6 
sin 6e~ %u) cos 9 



In Eq.(8) and hereafter, a quantity with subscript always denotes its value in vacuum. 
The mass matrix in the weak basis can be obtained from the mass matrix in the mass basis 
as 

MZ = U<,M?U^( mu mi2 V (9) 
y m 2 i m 22 / 

where M™ = &ia.g(m\,m%). 

In matter, a quantity A(t) = 2E u V(t), which is the amount of the increase in the effective 
mass squared in matter for u e , is inserted in the 1-1 element of the mass matrix 

. mii + A mio \ 
M w = | (10) 

m 2 i m 2 2 / 



with A = 2E U V = 2\f2GpE v N e . (We ignore the contributions from neutral current interac- 
tions.) The Gf, E v , and N e represent, respectively, the Fermi coupling constant, the energy 
of the neutrino under consideration, and the electron number density in a medium. Let us 
denote the effective mixing matrix in matter as 

/ cos 9 sm6e io \ , s 

V — sm 9e %w cos 9 J 

This effective mixing matrix U e ff diagonalizes the effective mass matrix in the weak basis, i.e., 
U\ffM w U e ff = M m , where M m = diag(/i 2 , /^|). The effective masses ^1,^2 and the effective 
mixing angle 9 are well known and are given by (For recent reviews, see for example, [0,pi|) 



A = + ~y(A-At cos 29Y + A* sin 2 29, (12) 

Hi = + \ + \ " A 2 , cos 29? + At sin 2 29, (13) 

^ 29 = (14) 

1 -A/A 2 , cos 20' v ; 

where Ag = m 2 — m 2 . In the two generation case, a simple algebra shows that the effective 
phase Q remains unaffected by matter, i.e., u = u. However, this is not the case for the three 
generation case as well as for the two generation case of the dynamical mixixing matrix, as 
we will see later. The above information is sufficient to describe the adiabatic conversion 
processes of u e in the Sun. 



B. Three Generation Case 

Now let us consider the case of three generations of neutrinos by generalizing the case 
discussed above. The representation of the unitary mixing matrix |T9|-pT| depends on the 
neutrino type. Majorana neutrinos can have three non-zero CP phases in vacuum whereas 
Dirac neutrinos can have only one CP phase. For simplicity, we will consider the case where 
only one non-zero CP phase appears in vacuum. The result which we will present below can 
easily be generalized to the case with more than one vacuum CP phases. 

We have found it convenient to describe the neutrino mixing with the use of the modified 
Maiani representation as advocated in Particle Data Group E3] given by 
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( C12C13 S12C13 Si3e lS \ 

— S 12C23 — c 12-5l3S23e J<S C12C23 — Sl2'Sl3' S 23 e4<5 C13S23 
Y S12S23 - C 12 S 13 C 23 e lS -C12S23 - Si2Si3C23e J ' 5 C13C23 J 



(15) 



where = siaOij, Cy = cos% . (The mixing angles % and phase 5 in Eq. (|T5| ) are all 
vacuum values although they do not have the subscript 0.) 

The mass matrix in the weak basis is related to the mass matrix in the mass basis as 

' mn run mi3 ' 

77721 m 2 2 77723 > 

V rn 31 m 32 777 33 J 



M™ = U M™Ul 



where M™ = diag(m 2 , m|, 7773). In a way analogous to the two generation case, the matter 
effect is inserted in the mass matrix: 

( 777ii + A 777i2 "^13 ^ 
M W = m 2 l m 2 2 77723 • (17) 

y m 3 i m 32 m 33 y 

Now, solving the problem in matter amounts to diagonalizing the Hamiltonian or solving 
the time evolution matrix in Eq. (jl]). In the adiabatic case or in the case of a uniform 
matter density, it is equivalent to finding the effective mass matrix M m and the effective 
mixing matrix U e ff from the equation U\jfM w U e ff = M m , where M m = diag(/xf, n% /•*!)■ 
The eigenvalues are 

r p + 2n. b 



i4 



2a 1/3 cos i 
2a 1/3 cos i 
2a 1/s cos i 



where 



a = \jq l + \p 6 + q z 



3 3 : 
(3-2n b 

3 ' ~ 3' 

r 3' 



(3 = arctan i 



(18) 
(19) 
(20) 



|p3 _|_ g2 
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9 : 



6c b 3 
~2 + 6~ ~ 27' 



with 6 = — TrM™, c = ^(TrM" 1 ) 2 - Tr(M w ) 2 ], and d = -detikT. 

Let us denote the eigenvector corresponding to the third eigenvalue as a column vector 
(ci,c 2 ,c 3 ) T , where T denotes transpose. Solving the equations for Ci,c 2 and c 3 , we obtain 
the ratios 

C2 _ (m 33 - fiQ 77721 ~ m 2 3in 31 

Ci 7772 3 777 3 2 - (m 22 ~ £*!) (t77 33 - [1%) ' 



(21) 



C3 = (m 22 - fiQm 31 - m 32 m 21 

ci m 23 m 32 - (m 22 - (m 33 - //§) ' 

In matter, the phases of c 2 /ci and c 3 /ci are no longer the same due to the presence of 
the non-zero vacuum CP phase 5. This is the reason why we are forced to introduce two 
additional phases for the mixing matrix to be unitary and self-consistent. Thus, we can 
not use the same form as the modified Maiani representation given in Eq. (]X5f) in order to 
describe the new effective matrix in matter. 

Let us take a representation in matter which has two additional induced phases, £ and r\ 
which appear with 9 23 and #i 2 , respectively, 



U, 



eff 



( c 12 c l3 s 12 c 13 e lv s 13 e~ lS \ 

-s 12 c 23 e~ lri - c 12 s 13 s 23 e liS -^ c 12 c 23 - s u s 13 s 23 e l{E+r >-0 c 13 s 23 e~^ 
V s 12 s 23 e~ 1 ^-^ - ci 2 si 3 c 23 e^ -c 12 s 23 e^ - s 12 s 13 c 23 e l( - s+v) c 13 c 23 ) 



(23) 



where = sin^, c\j = cos9y. A barred quantity denotes its effective value in matter. 
The two additional phases, £ and rj, are induced in matter due to the presence of the original 
non-zero CP phase in vacuum. With the representation Eq. (|23f) , we can obtain all the 
effective values in matter. 

The two effective mixing phases are given by 

tan*= jr! C3 M , ( 24 ) 



Re(c 3 /ci) 

M J_ O = £j2Z£0 (25) 
Re(c 2 /ci) 



and the effective mixing angles are 



tang 13 = (26) 
\J\c 2 /ci\ 2 + |c 3 /ci| 2 

COS 023 = , (27) 
\j\c 2 /ci\ 2 + I c 3 /ci 1 2 

In order to obtain the remaining mixing angle 9\ 2 and the phase rj, we must solve another 
eigenvector for an eigenvalue, for example, ji\. Let us denote the eigenvector as a column 
vector (di, d 2 , d 3 ) T . Solving the equations for di, d 2 and d 3 , we obtain the ratio of d 2 to d\ 

<k _ (m 33 - ji\)m 21 - m 23 m 31 

di m 23 m 32 - (m 22 - ji\) (m 33 - /if) ' 

After comparing with the ratios of the 1-1 and 1-2 mixing matrix elements, we obtain the 
effective phase rj and the effective mixing angle 6\ 2 as 



tan _ tan gig sin fl 23 sin(£ - Q + Im(d 2 /di) ^ 
tan 13 sin # 23 cos(<5 — £) + Re(c? 2 /^i) ' 

tan^i2 = — (tan 0i3 sin 2 3 cos(5 — + Re(d 2 /di)) = — — — . (30) 

cos 6^3 COS T] 



From Eqs. (^J) and (|29D , we can see that the two additional induced phases £ and rj are 



generated by the non-zero vacuum CP phase 5. That is, if we set 5 = 0, then Eq. fl24|) 
shows that the effective phase 5 also vanishes because c 3 /ci cannot have an imaginary part. 
Similarly, it is easy to see that the two additional phases £ and rj vanish if 5 = 0. We 
emphasize that the two phases £ and rj are not new and independent phases, but instead 
they are related to the non-zero vacuum CP phase 5. 

In order to exhibit this property quantitatively, we use the solar neutrinos as an example. 
Most solar neutrinos are produced in the region x < 0.25 where x is the ratio of the distance 
R from the center of the Sun and the radius of the Sun R & . In the Sun, the electron number 



density distribution ppfl is well fitted by \og w (N e (x)/N A ) = 2.32-4.17x-0.000125/[x 2 + 



0.5 2 ] where Na is Avogadro's number per cm 3 . The quantity A in the center of the Sun is 
A = 3.17 x 10" 5 (E l/ /MeV)eV 2 . With the solar neutrino energy E v < 14MeV, the largest 
value of A in the Sun is of the order of 10 _3 eV 2 . Therefore, if we assume small vacuum 
mixing angles, there are two resonance points (0y = vr/4) inside the Sun when the mass 
squared difference A 31 = m 2 — m 2 is less than 10 _3 eV 2 . 

Fig. 1 and Fig. 2 show the typical behavior of the three effective mixing angles and three 
phases in the case where the neutrino mass squared differences are taken as A 2 i ~ 10 _4 eV 2 , 
A 31 ~ 10 _2 eV 2 . In Fig. 1, we have taken the vacuum mixing angles #i 2 = #i 3 = # 23 = 0.05 
and the vacuum CP phase 5 = 0.5. Therefore, we can see that, in contrast to the case 
of the effective mixing angles and effective masses which are dramatically affected by the 
medium, the CP phase 5 remains unchanged, and £ and rj are induced for large values of 
A in the adiabatic processes. Fig. 2 is for the vacuum mixing angles Qyi = # 23 = 7r/4, 
#13 = 0.6155, and the vacuum CP phase 5 = n/2. This is the maximal mixing case. In Fig. 
2, the effective mixing angles 9 12 and 9 13 change from 9 12 = 7r/4 and #13 = 0.6155 in vacuum 
to 9 12 ~ 9 13 ~ 7r/2 in the solar interior, whereas the mixing angle 6? 2 3 and the phase 5 also 
practically remain unchanged. Again, the induced £ and rj appear when A becomes large. 
This confirms the assertion made in P5] , P6]| that the matter effects are practically absent in 
the maximal mixing case. 

In general, the non-zero CP phase contributes to the transition probability among the 
flavor neutrino states. For simplicity, we consider the averaged survival probability of u e in 
the adiabatic approximation. The survival probability is given by 



< P(y e — > v e \ adiabatic) > 



(C12C13C12C13) 2 + (S12C13S12C13) 2 + O13S13) 2 



+2C12C13S12C13C12C13S12C13 COS Tj 

+2C12C13S13C12C13S13 cos(5 - 6) 

+2S12C13S13S12C13S13 cos(?7 + 6-6). (31) 



Although three phases are present in Eq. (BTf), since they appear as cos^,cos(5 — r]), and 



cos (rj + 6 — 6), all of which are close to unity, the CP phases are almost impossible to detect 
in practice. (Note in the figure that 6 ~ 6 and £, 77 ~ 0). 

IV. DYNAMICAL NEUTRINO MIXING MATRICES 

A. Two Generation Case 

Now let us study the dynamical mixing matrix in detail. We will first consider the two 
generation case. From the definition of the dynamical mixing matrix given by Eq. ([5]), we 
can write the differential matrix equation in a 2 x 2 matrix form. As it becomes clear later, 
at least one non-zero phase is necessary independently of the neutrino type even in the two 
generation case. With the following notations 

(cos?/' sin^e"* \ 
— sm-^e cos ip I 

M e =\~ L " I , (33) 



and 



--r^ 2 + / (t) ! sin2 > (*> 

sin 20 cos 20 

Eq. (Q) leads to the two coupled equations for ip and a 

A 2 

ib = ^ sin 20 sin a, (35) 

A 2 n A 2 cos2#-A(t) , s 

a = sin 20 cos a cot 2^ + — — — . (36) 

2E 2E 

Here, it is easy to see that if there is no phase (i.e., a = 0), it is impossible to describe the 
case where A changes in time since with a = Eq. (|3~6"D is valid only at a fixed value of 



t. In the above, the mass matrix in the weak basis in Eq. ([34]) has been obtained from the 
mass matrix Eq. (|l0| ) by subtracting (m\ + rn^)/2 from the diagonal elements in Eq. (fL0|). 
In addition, we obtain the energy eigenvalues, £\{t) and £T 2 (i), 

A 2 A 2 
£ l (f) = — ^ cos 26 + A(t) - sin 26 cos ata.ni/j, (37) 

A 2 A 2 

£ 2 {t) = cos 26 - -y sin 26 cos a tan ip. (38) 

Equations ( p5|) and (|36|) are coupled non-linear differential equations. Similar equations 
have previously been obtained in [fj7| as a means to describe the non-adiabatic process 



without the use of the Landau-Zener transition probability. The initial conditions for ip and 
a are fixed by the vacuum values. Given A(t), we can solve the equations numerically. For 
E/(m,2 — mf) — > oo, ip(t) and a(t) approach the vacuum values, which can be seen from 
the usual treatment with the Landau-Zener formula. Fig. 3 (a) (Fig. 3 (b)) shows the 
dynamical mixing angle and phase in the case of the neutrino energy E u = 10 MeV, the 
vacuum mixing ang le 6 = 0.1, and A 2 = 10~ 4 eV 2 (A 2 = 10" 6 eV 2 ). The difference between 
the effective and dynamical mixing angles depends sensitively on Aq. Fig. 3 (c) shows the 
dynamical mixing angle and phase in the maximal mixing case of two neutrinos with the 
neutrino energy E u = 10 MeV and the mass squared difference Aq = 10 _6 eV 2 . The abscissa 
x indicates the ratio of the distance of interest from the surface of the Sun to the solar 
radius, i.e., x = corresponds to the surface of the Sun and x = 1 the center of the Sun. 
The sign of the dynamical phase a depends on the direction of the neutrino motion, i.e., 
the sign of a is reversed as v e travels inwards from the solar surface. The behavior of a as 
a function of x is quite remarkable as x — >• 1. 

The time averaged survival probability for v e is given by 

< P(u e -> u e ) >= cos 2 i/j cos 2 6 + sin 2 ip sin 2 6. (39) 



We emphasize here that although the expression in Eq. (|39"D is identical to the usual survival 
probability for the adiabatic process when ip is replaced by the effective mixing angle 6, it is 
an exact one which can describe the non-adiabatic process without using the Landau-Zener 
transition probability. In Eq. fl5§D the dynamical phase a has disappeared in the process 
of time average. Fig. 4 shows the survival probabilities given by Eq. (^) as functions of 
E v /Aq for three examples with the vacuum mixing angles 6 = 0.01, 6 = 0.1 and 6 = 7r/4 
(maximal mixing in two generations). The transition probability for the maximal mixing 
case is practically constant at 0.5, in the region of 10~ 10 < E u /Aq < 10~ 15 and the neutrino 
oscillation in matter cannot be differentiated from the pure vacuum oscillation. However it 



does not mean that the dynamical mixing angle ip remains constant in matter. The value of 
the angle ip approaches to 71/ 2 near the center of the Sun. In Fig. 4 the survival probabilities 
calculated using the Landau-Zener transition probability (dotted lines for 9 = 0.01, 9 = 0.1) 
are compared with the exact results based on the use of dynamical mixing angle. 



B. Three Generation Case 



In the case of three generations, there are six coupled differential equations and three 
relations to be solved. The unknown variables are three dynamical mixing angles, three 
dynamical phases and three energy eigenvalues. Let us write the dynamical neutrino mixing 
matrix in the three generation as 



U, 



dyn 



(l 







C23 



\ 



\0 -s 23 e i7 c 23 / 



V 



C13 s 13 e-tf \ ( 
1 

-s 13 e i(3 C13 



J 



C12 s X2 e 

ia 



\ 



s 12 e 




C12 








(40) 



where c^- 



cos ijjij,Sij = sinipij with a, (3 and 7 denoting the dynamical CP phases in 
matter. From the off-diagonal elements of Eq. (|5]) we can obtain the following differential 
equations for three mixing angles and three mixing phases, 



B 





(i-n) 






fk.3 
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2E 
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{ 7 ) 




x 6 matrix B is given by 


-s a 


C?2 s 13-S/3- 7 — sf 2 


c a 


C 12 5 


13 c /3-7 + sf 2 


C 12 S p 


S12C13SQ 



12 
12 

13 

Im(f4 n M»fV) 

13 

Re([/J yn M-[/ ( 



Im(E/J_AW, 



dyn ) 23 
dyn) 23 / 



"7 
-7 







C12C/3 



Cl2Cl3S 7 Si2Sq, + /3 
V Ci2Ci3C 7 Si 2 C a+ ^ 



— Si 2 Ci3C a _ 7 






-Sl2Cl2C a 
-Si2Ci 2 S Q 








-s 12 c 12 sj 3 c a 




2 

— s 12Cl2-S 13 S Q , 


b 26 


C12S13C13C/3 


h& 


-C12S13C13S/3 


bid 


•5l2Sl3Ci3CQ, +i a 


he 


— Sl2Sl3 c 13 s aH-/3 


bm j 



(41) 



(42) 



where 



blB = Sl2Cl 2 (l + S 2 13 )sl 3 C a + sl 2 S 13 S23C 2 3C2a+(3-y ~ C 2 12 S 13 S 23 C 23 Cf3- 

b 2 6 — Si 2 C\ 2 (l + sl 3 )sl 3 s a + s 2 2 Si 3 s 23 c 23 s 2a+ p^^ + C^ 2 S\ 3 S 23 C 23 Sp. 



■11 



-7 J 



he = -Sl 2 C 13 S 23 C 23 C a - y - C 12 S 13 C 13 sl 3 Cp, 
&46 = -Sl2Cl 3 S23C23S Q - 7 + C 12 S 13 C 13 sl 3 Sp 



1 ! 



&56 — — ■Si2Si3Ci3S 23 C a+( g + C12C13S23C23C, 
^66 = s 12 s 13 c 13 s 23 s cH-/3 — Ci2Ci3S23C 2 3S 7 . 

The energy eigenvalues are 

= (I^A^E^ii - 2E{as\ 2 - (3c 2 12 s 2 13 - (43) 

7[( s 12 - C 12 S 13) S 23 _ 2si2Ci2Si 3 S 2 3C23Cq+/3- 7 ] - 2V>23Sl2Ci2Si 3 S a+/3 - 7 }, 

= (U\ yn M w U dyn ) 22 + 2E{as\ 2 + (3s 2 12 s 2 l3 + 

7[( C 12 — S 12 S 13) S 23 ~ s 12Cl2'Sl3' S 23 c 23 c a+/3- 7 ] — V , 23 s 12Cl2-5i3S Q , +( g_ 7 }, (44) 

= {U\ yn M w U dyn ) 33 - 2E0s\ 3 + ic\ 3 s 2 23 }. (45) 
The averaged transition probability is of the following form: 

< P{v e -> i/ e ) >= £ K^VM^okl 2 - (46) 

i=l 

In the three generation case, the dynamical phases contribute through (U^ yn )a to the transi- 
tion probability in contrast to the two generation case. The formulas in the three generation 
case are reduced to the corresponding expressions in the two generation case, when we take 
#13 = #23 = and m 3 i = = 0(i = 1,2,3). In the case of three generation of neutri- 
nos, however, it is prohibitively difficult to solve the coupled differential equation (|4l|) - ( f45"D 
numerically and such a study is beyond the scope of the present work. 



V. CONCLUSION 



In this paper we have introduced the effective mixing matrix that relates weak eigenstates 
to mass eigenstates, and the dynamical mixing matrix that connects weak eigenstates and 
energy eigenstates. Using these definitions, we have studied the effective mixing matrices 
for the two and three flavor mixing cases, with emphasis on the nature and behavior of 
the CP phases in matter. It has been shown that in the case of three generation mixing, 
two additional induced phases appear when a vacuum CP phase is non-zero. Analysis of the 
dynamical mixing matrices for the two and three flavor mixing cases has also been presented. 



In contrast to the effective mixing matrix of the two flavor mixing, in the dynamical mixing 
matrix, it is necessary to introduce one additional phase in a medium in order to be self 
consistent. This does not depend on the neutrino type. We have used the solar neutrinos 
to demonstrate the differences involved in the behavior of the CP phases in matter, and the 
differences between the effective and dynamical mixing matrices. 
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Figure Captions 



Fig. 1 The effective masses, mixing angles, and CP phases in matter as functions of A with 
m \ = l(r 6 eV 2 ,m 2 = 10~ 4 eV 2 ,m 2 = 10" 2 eV 2 , mixing angles 9 12 = 9 13 = 9 23 = 0.05 
and vacuum CP phase 5 = 0.5. 

Fig. 2 The effective masses, mixing angles, and CP phases in matter as functions of A 
for the maximal mixing case: m\ = 10~ 6 eV 2 ,m 2 = 10 _4 eV 2 ,m| = 10~ 2 eV 2 ,vacuum 
mixing angle 9 12 = 7r/4, 9 13 = 0.6155, 9 23 = 7r/4, and vacuum CP phase 5 = ir/2. 

Fig. 3 The effective (9) and dynamical (ip) mixing angles and the dynamical mixing phase, 
a in the case of two generations. The parameters used are E v = 10 MeV with (a) 
the vacuum mixing angle 9 = 0.1, Aq = 10~ 4 eV 2 , (b) the vacuum mixing angle 
9 = 0.1, Aq = 10~ 6 eV 2 , and (c) the maximal mixing with two generations 9 = ir/2, 
Aq = 10 _6 eV 2 . The abscissa x indicates the ratio of distance from the solar surface 
to the radius of the Sun. The value x = means the surface of the Sun and x — 1 
the center of the Sun. The sign of the phase a is reversed when the neutrinos travel 
inwards. 

Fig. 4 The u e survival probabilities in the two generation case as functions of E„/Al for 
the vacuum mixing angle 9 = 0.01, 9 = 0.1, and 9 = 7r/4 (maximal mixing). The 
solid curves are the results with the use of the dynamical mixing angles whereas the 
dotted curves (9 = 0.01, 9 = 0.1) are those with the use of the standard two generation 
formulas with the Landau-Zener transition probabilities. The Landau-Zener formula 
cannot be applied to the case of large mixing, specifically to the maximal mixing 
mixing case. 
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